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Abstract 

Motivated by recent measurements which strongly support a non-zero reactor mixing angle ^13, we study a 
deviation from 6*3 neutrino discrete symmetry by explicitly breaking the neutrino mass matrix with a general 
retrocirculant matrix. We show that non-zero ^13 and non-zero CP violation parameter Jcp arise due to the 
difference between 2/2 and 2/3. We demonstrate that it is possible to obtain the experimentally favored results for 
neutrino masses and mixing angles from this mass matrix. Furthermore, we estimate the effective masses mp and 
mpp and total neutrino mass l"^'! predicted by this mass matrix. 
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Experiments using solar, atmospheric and reactor neutrinos have made considerable progress in establishing 
two different mass squared differences ( Am^j^ and Amgj^ ) and two large mixing angles ( 9i2 and 623 ) in the 
lepton sector. Recently MINOS [Il[2], T2K 0, Double CHOOZ Il],DayaBay [S| and RENO ^5" have revealed 
that the reactor mixing angle ^13 is not only nonzero but relatively large. 

The phenomenon of neutrino mixing can be simply described by the Pontecorvo-Maki-Nakagawa-Sakata ( 
PMNS ) neutrino mixing matrix Vpmns 13 j which links the neutrino flavor cigenstatcs i/e,i';^,t^r to the mass 
eigenstates i^i , , ^^2 ■ 
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In the standard parametrization used by the Particle Data Group ( PDG ), the PMSN matrix is expressed by 
three mixing angles 9i2, 023 and ^13 and one intrinsic CP violating phase i5 for Dirac neutrinos, 

(C12C13 S12C13 sise'"-^ \ 

-S12C23 - Ci2S23Sl3e"^ C12C23 - Si2S23Sl3e"^ S23C13 -PMaj (2) 
S12S23 - Ci2C23Sl3e"^ -C12S23 - 5120235136"^ C23C13 / 

where Cij = cos^y, Sij = sin^^j and P^aj is a diagonal matrix with Majorana CP violating phases. 
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The three mixing angles are related to the moduli of the elements of the PMNS mixing matrix as : 

sin' = 1 143 1' 
sm Oi2 — z 
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The well-known Tri-Bimaximal (TBM) mixing pattern, which corresponds to ^13 = 0, 623 = ±^ and 
6*12 = sin^^ ("Ta)' attracted a degree of attention in the literature because it suggests some underlying 
flavor symmetry among lepton's generations. This flavor symmetry is expected to explain the mass spectrum and 
neutrino mixing pattern. 



The TBM form is 
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Models based on discrete symmetries were successful in reproducing this matrix. Among a number of inter- 
esting discrete flavor symmetries discussed in the literature, the S3 symmetry which is the permutation group of 
three objects, is the simplest [S]. ^3 is the smallest non-Abelian discrete group. 
The three dimensional reducible representations of all 15*3 group elements are : 



5(1) =010 ; 5(12) ^ 100 ; S^''^ 






^(23) =001 ; ^(123) ^ 1 I ; ^("2) = ( 1 | . (5) 
\010/ \010y 

The most general neutrino mass matrix invariant under 6*3 is : 

MO = a5(i)+/3(5(i2)+5(i3)+5(23)j (6) 

where a and /3 are, in general, complex numbers. 

In the basis where the charged lepton mass matrix is diagonal, the TBM mixing matrix diagonalizes the 
neutrino matrix M^. 

fa \ 

Vo'^Myo = a + /3 . (7) 

\ a J 

This matrix leads to two degenerate masses, namely mi and 7713. However, this is not correct experimentally. To 
overcome this problem, it was suggested in [3] that in fact the three masses are degenerate by letting the complex 
number a lies in the third quadrant. 



-i\a\e~'T^ for < < tt (8) 
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and taking /3 as real number such that : 



P = • (9) 



In this work, wc consider the neutrino matrix M^, which is invariant under 6*3, as zeroth order with degenerate 
masses and TBM mixing angles. Non-degenerate mass spectrum and non-zero ^13 were realized in [9]-|10j by 
introducing small perturbations that violate ^3 symmetry. 

Here, we investigate the phenomenological consequences of the deviation from an exact symmetry by 
explicitly breaking the neutrino mass matrix with a general retrocirculant matrix : 

/ yi y2 Vs \ 

AAU = -a 2/2 ys yi (10) 
\ 2/3 yi y2 j 

where the dimensionless parameters yi = |yi|e"^'/^ are complex numbers with magnitude less than one and 
< 5: We also consider the charged lepton to be diagonal so the leptonic mixing solely comes from the 
neutrino sector. It is easy to see that AMj, can be written as a linear combination of S'(^'^\ 5^^^^ and S'^^^^ as : 

As a result the broken neutrino matrix M^, becomes : 

a + P - ayi jS ~ a y2 P - a yz 

13 - ay2 a + (3- ayz P - a yi | . (12) 
P - ay:i P - a yi a + P ~ a y2 

The eigenvalues of the above matrix are : 



mi = a~a yl + yl + yl- yiy2 - 2/12/3 - 2/22/3 
1712 = a + 3/? - a (j/i -I- 2/2 + 2/3) 



m3 = a + a Y 2/1 + 2/2 + yi ~ J/iya - 2/i2/3 - 2/22/3 • (13) 

The matrix M^, is called magic mass matrix since every row and column add up to the same value which is 
1712 for this case. It implies that this mass matrix has a trimaximal eigenvector ( -1= , -1= , --i= j [IT]. Moreover 
the fi — T symmetry corresponding to ^13 = and 623 ~ ±j is broken for M^, 

[M.,5(23)] = _«(y,_y3) (5(123)_^(132)^ (14) 

due to the difference between 2/2 and 2/3. 

Interestingly, by rotating the above matrix by the TBM mixing matrix Vq, we get : 

V - f (22/1 - 2/2 - 2/3) ^ (2/2 - 2/3) 

Vo^M.Vo - I a + 3/3 - a (2/1 + 2/2 + 2/3) | . (15) 

^(2/2-2/3) a H- f (22/1 - 2/2 - 2/3) 
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As a consequence, the neutrino matrix M^, is diagonalized by the total unitary matrix V = VoUPMaj- The 
mixing matrix U is given by : 

cos 9 e"^ sin ( 
U = \ 1 ) . (16) 

-e-*'' sin 61 008 6* 



A straightforward calculation yields that the angle 9 and the CP-phase S are 

tan 26* 
tan 5 
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where 
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The explicit expression of V is 
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(19) 



We immediately obtain the mixing angles : 



sin^ 6*13 



sm 9i2 
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Z + 2VX^TW+Z^ 



(20) 



which for ly^l = |y2| = give the TBM mixing angles. 



The strength of CP violation in neutrino oscillations is described by the Jarlskog rephasing invariant parameter. 
It is given by : 

1 

>2) - 



JCP = /™(K2^p3K3^; 
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It proves convenient to use the ratio i?^ 



Am 



r, to write the mass squared differences as : 



Amji = 2|api?^V^2 + y2 + Z2 
Am^i = 2|a|V^2 + y2 + ^2 



'32 



2|ap(l - R,) + y2 + ^2 



(21) 



(22) 
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Ill order to confront the above neutrino matrix with the experimental observations, we use the constraints on 
neutrino parameters at la and Scr |12| . 



Parameter 


best fit 


2cr 


3(T 


Am^JlO-^'ey-'] 


7.62 


7.27-8.01 


7.12-8.20 


Am^JlO-^ey'-^] 


2.55 


2.38-2.68 


2.31 - 2.74 




-2.43 


-(2.29 - 2.58) 


-(2.21 - 2.64) 


sin"' 6*12 


0.320 


0.29-0.35 


0.27-0.37 


sin^ 6*23 


0.613 


0.38-0.66 


0.36-0.68 




0.600 


0.39-0.65 


0.37-0.67 


sin'' 6*13 


0.0246 


0.019-0.030 


0.017-0.033 




0.0250 


0.020 - 0.030 


0.017-0.033 




O.SOtt 


0- 27r 


- 27r 




-0.037r 


0-27r 


- 27r 



Table 1: Global oscillation analysis with best fit for Am^^, Amgj^, sin^ 6*12, sin^ 6*23, sin^ 6*13 and 6 the upper and/or 
lower corresponds to normal and/or inverted neutrino mass hierarchy. 



For numerical analysis, we use the 3a ranges of neutrino oscillation parameters. The neutrino mass matrix 
depends on |a|, jj/ij, |j/2|, lual and the phases 'ip,ipi,ip2 and ip^. 

Since there are many unknown parameters, we consider a particular set of those parameters and show how the 
measured values of neutrino experiments can be accommodated in our neutrino mass matrix M^. For simplicity, 
we take ip2 = ~ 0, as inputs. 

To see how the non-zero value of the phase ipi can lift the degeneracy between mi and m2, we notice that for 
2/3 = 2/2=0, the mass squared differences become : 

A 2 /II |2| I ■ V' . -ip-fl 

Am^i = 4|a| |yi|sin-sin 

Am^i = 4|a|2|yi|cos^ 

Aml^ = 4|anyi|cos|cos^^^ . (23) 

It is worthwhile to remark that for ipi = ip, the masses mi and m2 are degenerate and |m3| > |m2|. Now to 
separately obtain a non-zero mass squared difference Am^i which is smaller than Am|i, we introduce a small 
phase difference ei = ~ (pi between the phases tp and (pi . Such a small phase difference will be responsible for 
lifting the mass degeneracy between the first and second generation. 

To sec the behavior of the mass eigenvalues with respect to ip and other observables, we take ei = 6° as a 
typical value with \a\ = 0.1 eV, \yi\ = 5 x 10~^, \y2\ ~ 1.3 x 10~^ and |j/3| = 1.2 x 10~^. Figure 1 shows the 
variation of the masses \mi\ as a function of the phase ip- One clearly observes that Fig.l suggests a normal 
hierarchical ordering pattern for ip > 70°. 

Based on the expression of the neutrino mass squared differences, we numerically scan over a broader range 
of \yi\ (0.05 < \yi\ < 0.4) to obtain the restriction on the parameter space of \y2\ and ly^]. We have plotted in 
Fig. 2 the contour plots of the mass squared differences in the two dimensional parameter spaces (|2/2|j I2/3I) where 
we take 4' — 120°, ei = 6° and |a| = 0.1 eV as inputs. 



5 



0.115 r 




Figure 1: Neutrino's masses versus ijj. 



From Fig. 2, we obtain the allowed ranges of \y2\ and jj/aj for the normal mass hierarchy, 

|y,|<0.02 fori = 2, 3 



(24) 




The ratio R^, for normal mass hierarchy ordering in the 3cr allowed range is : 

= (2.99t°j2) X 10-2 . (25) 

Fig. 3 shows contour plot of the ratio R^, in the two parameter space (ip, ei = ip ~ tpi) where \yi\ = 0.05, \y2\ = 
1.3 X 10^'^ and jyaj = 1.2 x 10^'^. As expected, the ratio Ri, at 3(t indicates that ip has to be greater than 70°. 



The departure of the mixing angles from TBM mixing angles depend on the phase ei and the two parameters 
|?/2| and lyaj. Scanning over \y2\ and jj/sl within their allowed ranges ( |y2|, jj/sl < 0.02 ), we investigate how 
a non-zero 813 can be obtained for normal mass hierarchy. As a result of numerical analysis, contour plots in 
the (-0, \yi\) parameter plane of the mixing angles ^13, 612 and 623 are shown in Fig. 4 using the experimental 
constraints on the measured angles. 

The CP violation parameter Jcp, which is directly related to the Dirac phase 6, arises due to non-zero value 
of the difference {ysl — \y2\ in the neutrino matrix M^. We have plotted in Fig. 5, Jcp with respect to using the 
allowed region of \yi. It leads to values of \ Jc'p\ around 10"'^. Non-vanishing \ Jcp\ will be explored by the next 
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Figure 3: Contour plot of the ratio in the parameter space (i/;, ei)- 




Figure 4: Contour plots of reactor angle 6*13 (left panel), solar angle 6*12 (middle panel ) and atmospheric angle 
023 ( right panel ) in the parameter space |yi| and '(/'• 




Figure 5: CP-violating parameter Jcp versus "0. 
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generation high performance long-basehne neutrino experiments. 



The absolute neutrino mass scale can be probed by non-oscillatory neutrino experiments. Cosmology is sensi- 
tive to the sum of neutrino masses l"^*!- The beta decay endpoint measurements probe the so-called effective 
electron neutrino mass m^. The rate of the neutrinolcss double beta decay depends on the effective Majorana 
mass of the electron neutrino m^^. 

Both mpp and and the sum of neutrino masses are given by : 

^fSP = I H ^i^ei I 



'^0 = jXI'™?!^^^! 



m\ = -I- |m2| -f Imsl . (26) 



From the neutrino mass M^,, the effective Majorana masses m/jp and mp can be written as 



mpfs = ^ W5 + 4cosV'-6|yi| ( sin "^^ +2sin^ ) +9\yi\ 




Present cosmological constraints on the sum of neutrino masses X]!"'^^ ^^'^ ^1^*^ range 0.44 — 0.76 eV [13) . 
The Mainz |14j and Troitsk |15j experiments on the high precision measurement of the end-point part of the 
/?-spectrum of decay found the 95% C.L. upper bounds mp < 2.3 eV ( Mainz ) and mp < 2.1 eV ( Troitsk). 
Experimental bound on mpp is below 0.36 eV [16j . 




Figure 6: mfj^mpp and ^ \mi\ versus ip. 

Fig. 6 gives the effective electron neutrino mass m^, the effective Majorana mass mpp and the sum of neutrino 
masses ^ \mi\ with respect to ip for ei =6°. It shows the predicted m^, mpp and ^ \mi\ are well below the 
experimental bounds. The magnitude of mpp increases with larger \yi\ values. 
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